

“5 
<3 * 


-Si ^ 

.*= o> 


S S3 

fci 

«3 a? 
« 1= 

-o «? 


s 3 

c: at 3 
=1 ^ -■ 

811 

sn 

o o S 


© “ 3 

s e ^ 

2 S C 

g g ; 

* £ £ 


S-193 IMPULSE RESPONSE CROSS CORRELATION 
(Report Covering Sept. '73-March '74) 

by 

Dr. H. S. Hayre 


WORK DONE UNDER 
NASA Contract #NAS 9-13462 
March 29, 1974 
TR— 74-11 



s1 H-28885 


Unclas 

CSCL 05B g3 ^ 1 3 00652 



TECHNICAL MONITOR 
T.T. WHITE MAIL CODE TF6 
P.i. MANAGEMENT OFFICE 
NASA LBJ SPACE CENTER 
HOUSTON, TEXAS 77058 


UNIVERSITY OF HOUSTON 

AVE PROPGAT I ON/LAB. -ELECTRICAL ENGINEERING DEPT. 
CULLEN .COLLEGE OF ENGINEERING 
/HOU ST ON, TEXAS 77004 

/Tel.s -4c 713-749-4503/4535/4412 


i 


Abstract / Summary 

The impulse response of terrain to Altimeter (S 193) 

Model I, periodic pulse train signal and its verification consists 
of ground truth and terrain topographic data collection as well 
as sky lab data procurement and computer analysis, and finally 
the cross correlation between the sky lab data derived impulse 
response and that obtained from the corresponding topographic 
data. The ground truth and topographic (elevation) data 
were collected and reduced for Oregon Coast (SL #2) and Colorado 
overpass (SL #3). Seven track sky lab data was replaced by nine 
track data by December, 1973 in order to save excessive computer 
;cost Retrieval of data usable in this study was initiated 
by late February 1974. Computer program s : for calculating the 
terrain impulse response from sky lab data were written, dry 
run, rewritten and modified* They are now considered operational* 
The matter of frequency and time resolution in computer calculations 
is still receiving more study. 

A significant result for this period was the realization 
that the phase information, normally lost in envelope detec— 
tion and sampling, could have been preserved if the samples 
Were taken at the peaks of the carrier frequency at intervals 
equal to a twice or an integral multiple of the integral number 
of its period. Thus the Skylab S 193 Model I altimeter data 
does not contain phase information in accordance with the tech- 
nical information obtained by us from NASA and Lockheed technical 
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personnel. Another significant discovery was to learn that the 
eight return pulse samples don't belong to the same pulse but 
represent an average over many pulses. 
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INTRODUCTION 

The work for the period September, 1973 through March 31, 

1974 is covered in this quarterly report, since the skylab 
data received until November 30, 1973 was not considered sufficient 
to yield any significant practical results. The reason was 
that skylab data received until November 30, 1973 was in seven 
track-magnetic tape format, and the reading of that data and 
reformating it to meet our Univac nine track requirement would 
have been excessively expensive and not possible within the 
budget provided. The request for nine track tape was processed 
by NASA-LBJ Space Center, and the first nine track data tape 
on S 193 Altimeter Mode 1 was received in late December^, 1973 . 
(Future reports shall be submitted quarterly.) Then it was 
possible to dry run our computer programs on that data to verify 
basic needs as to resolution in time and frequency domain for 
the terrain impulse response calculation from this data. 

Various aspects of our work under this contract with their 
present status are listed below under separate subtitles. 

1 . SKYLAB- ALTIMETER DATA 

(A) S 193 Altimeter Mode-. 1: Since the start of and receipt 

of the nine track magnetic tape Altimeter (S 193) Mode 1 data 
in December., the following-data tapes have been received: 

1. SL-2 Project No. 8552C Nasa tape # 906153 (sl9313-101-2-l" 

73 -: 

Covering Oregon Coast Time: 150-20-35-41 

* 

.Thru: 150-20-46-45 
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2. SL-3 Project No. 8550 Nasa tape # ( S _93_i00-01-12-72A) 

907259 

Covering Colorado Time: 215-17-54-13 

Thru: 215-18-16-29 

3. SL-2 Project No. Nasa tape # 906154/s 

Covering the Gulf Coast and Montana Time: 160-15-03-30 

160-15-07-14 

160-15-11-35 

160-15-18-4 

4. SL-3 Project No. 8550 Nasa tape # S193-13-099-01-42-73-B 

Covering Mexico and western part of Texas 

Time: 258-16-23-06 

258-16-45-05 
258-16-23-9 
258-16-18-4 

Since we are unable to associate tape numbers three and four 
above with preplanned flight and sites for Altimeter Mode 1 
overflight, although efforts to salvage usable information 
£ rom tapes is continuing to be made. Data from Tapes numbers 
one and two was taken over the sites assigned namely # 851167 
and # 398295 respectively, and hence his data is being analysed. 

(B) S’ 191 and 190 and Associated Data 

Other associated data including transparencies of color film 
shots covering two sites have been received for SL-3 corresponding 
to dat tape numbers two and four above since January 1974. 

An exact coordinate match is underlay forvisual verification - 
of data. 
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(C) General and Specific Information 

The product S073-5 for S193 Mode I has-been requested for sites 
covered for this study, since only eight point pulse data is 
available, and the resolution enhancement is a major concern 
as discussed below. 

Ground Truth Data 

The ground truth data consists of such factors as moisture, 
rain fall, foliage, season, man-nade effects such as timber 
cutting, mining or farming, forest fires etc v type of forestation, 
crops or grass etc., and topographical data. The status of this 
work to date of the Oregon and Colorado sites is as follows: 

1. Site visit completed 

2. Photographic and visual information on landscape, 
foliage, crops, or forestation, man-made effects collected 

3. Topographic maps have been read, and data collected 

4. Ground track of skylab has been delineated in order 
to pinpoint any special features along it 

The S-193 Altimeter Mode I data over the Texas Hats Site 
apparently was not taken during all skylab passes, and therefore 
the ground truth along the overpass over the planned Hats site has 
not been compiled, although topographic maps for this planned 
overpass were procurred. 



Topographic Data and Analysis 

The resolution of topographic data is 80 feet for all quad- 
rangles in the Colorado overpass except for the quadrangles of 
Nucula, Silverton, Wolf Creek, Charaa Peak, and Brazos Peak; for 
which the 208 feet resolution data will be interpolated to yield 
points 80 feet apart, resulting in expanding these sectional data 
points by a factor of approximately 2*5 . Similarly, the resolution 
for the topographical data for the Oregon overpass data is 208 
feet. 

This data is ready for correlation analysis using a moving 
window, which is equivalent to a period .of the S-193 pulse , so 
that the pathlength on the ground illuminated by the altimeter pulse 
could be the same as the width of the window in our data processing. 
The only delay in initiating this processing is the receipt of de- 
tailed information on the S-193 transmitted pulse shape, since the 
point of return form the ground shall depend on cT/2 where c is the 
velocity of light or 3 x 10 s meters/sec and T is the effective pulse 
width, which is approximately 72 to 100 nanoseconds. 

Furthermore, it has just come to our attention that in S-193 
pulse shape experiment Mode I, the eight samples are so taken that 
theses 

a. do not belong to one return pulse 

b. are not consecutive in information sense, except in some 
"average" sense. 

c. can not be said to represent a single foot print of S-193 
beam on the ground 



The Skylab Altimeter S-193 Mode I pulse data was taken using 
a pulse train with a nominal pulse width of 100 nanoseconds and a 
repetition rate of 250 pulses per second. Furthermore, the return 
pulse was sampled eight times during each pulse duration, and the 
sample spacing was 25 nanoseconds. This would yield a spatial sampl- 

i 

ing rate, for the terrain under the skylab, of a sample every 84,5 feet, 
since the skylab ground speed is approximately 4 miles per second or 
21 feet per millisecond. Therefore the pulse period is the equivalent 
of 4 x 21 = 84 feet on the ground. 

An exact detail of sampling method or routine has been requested 
from NASA for a clear understanding of which ground samples, the aver- 
age received pulse represents. 

Transmitted and Received Pulse 

In Mode I pulse shape S-193 Altimeter experiment, due to no de- 
tailed information on the shape of the transmitted pulse x(t), it 
was assumed that it was an ideally square pulse of 13.98 ghz carrier. 
Later on inquiries on the rise and decay time of this pulse led us to 
find out that the pulse shape is indeed not spuare , but somewhat 
gaussian; but again the details of the clipped gaussian pulse are not, 
to this date, available from either NASA or NASA Contractors, such as 
Lockheed Electronics and Research Triangle Inc. in North Carolina; 
although these are expected to be received before the next quarterly 
report is issued. 

The importance of the shape of this transmitted pulse is demon- 
strated by the Fourier Transforms given below; 
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Case I : Rectangular Pulse , 

If the envelope of the transmitted pulse x(t) has a period T and 
pulse width a, then its Fourier Transform pair would be given below: 




c 



^ ‘v*. 
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Futher more the power spectral density P (f) -• would be given as: 
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Case II : Gaussian [approximated .by (sin x/x) form ] 

If the envelope of the transmitted pulse x(t) is approximated 
by an gaussian curve, as understood by telephonic inquiry as of April 
1/ 1974, then its Fourier transform pair shall be given as: 
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Note again that the phase spectrum is zero here as well. 


Case III: Truncated Gaussian Pulse 


The real case of S-193 Mode I pulse is closest to a truncated 
gaussian pulse. For instance, let us say over, one period, one may 
write x(t) = A exp (-K 2 t 2 ) -T/2< t <T/2 and its generalized 

fourier coefficient C shall be defined as: 


and 
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Upon substitution of limits in the value of Cn, one obtains: 

Cn = (2/Tk) (n/16) 1/2 exp(Wn/2k)2[erf (kT/2-jWn/2k)+erf (kT/2+jWn/2k) 

-2erf (jWn/2k)J 
since erf (-x) = erf (+x) 

It is important to note that the power spectral density of x (-t) 

. 2 

is iC n l / and its phase spectrum is not the same as in cases I and II 

Case IV : Truncated Gaussian Pulse with Fast Drop Off , 

In this case the Fourier Transform integral for C n in the 
previous case may be modified to extend its upper and lower limits 
to infinity since the pulse is already assumed to have dropped to 

negligible and essentially zero amplitude before reaching the 

✓ 

limits of its period. Thus we would obtain: 

C n » x(f>n = (2/T) :?exp(-k2t2) cos w t dt 

= (2/T) (II/4k2 ) exp (-w2/4k^) 

Note that in this approximation the power spectral density of x (+) 
in a given period is essentially the same as that obtained in 
Case III but its phase spectrum is zero as was the case in Case I 
and Case II. 

Importance of Power Spectrum and Phase Angle Spectrum in Impulse 
Response Calculations . 

The results obtained above will now be shown to have sig- 
nificanCe in this study, where the power and phase spectrum 
of x(-t) is needed to calculate the impulse response of the terrain 
radiated by S-193 pulse. Let the received pulse be y(+), and its 



power spectral density Py (w) , the phase angle spectrum Py(w ) , and 
the corresponding values for x (-fc) as P x (w) an< ^ tfbc-Cw) • Then the 
Fourier Transform of h (*t) , the terrain impulse response is given 
by the following relationships: - 

Y(jw> = X(j*>> 

Py CW) = P„ CW) 

<pj& - 9x (u,) * ft Cu ° 

“ . 

and therefore 

Ktb), ± [ p h fe «* P . j <*«■*> A*- 

— «o 

In this connection the Fast Fourier Transform technique 
was used to solve for h(t) given the eight sampled values of x(^). 

A typical program for this effort is attached and it is the result 
of extensive computer analysis of data samples of the type available 
from S-193, with basic intent to enhance resolution of the h (-fe) . 

Most early programs yielded a quick drop-off of the values 
of h £-fc) for t greater than zero. In fact the first value, for the 
data shown on attached pages A1-A6 , would be 176*51 and the next 
would be 0*539. A special technique of padding data with optimum 
number of zeros has now yielded at least 4 to 8 point resolution. 

In this effort eight S-193 return pulse sample values are extended 
to 256 points by adding 248 zeros at the end, and assuming that the 
phase angle information for the eight points is unknown and assumed 
zero as is the case of envelope detection. 

Other programs with phase angle imput, with eight points 
of x(«t) were also tried and the results were significantly interesting 
A close analysis of this part of the work reveals that this system 
would have been much more effective for such work and all other 
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surface effects if the pulse samples were taken at the peaks of 
the carrier signal and synchronized to the starting of the pulse 
at its positive zero crossing, and thus preserving its phase informa- 
tion. 

Conclusion : 

Recommendation for a new simple technique has been advanced 
to enhance the capability of S-193 Mode I altimeter pulse shape ex- 
periment* Futher work on the transmitted pulse shape specification 
and its effect on the impulse response of terrain shall be reported 
in the next quarterly. Futhermore cross correlation between 
different skylab data and the impulse response calculated shall be 
reported in the next report. 



•. , ■ • ; '• :* . . .. 

OR * US SKYLAB.MA iN, .MAIN ' . 

F OR S0ll-04/04/74-l3:47:34 -<2.3) . 

main program 

STORAGE USED: CODE(l) 000277; PATAU!) 010103; BLANK COMmON<2> 000000 


EXTERNAL REFERENCES (BLOCK. NAME V-y ■/ 

■0003 . FOUR1 ; 'O; : A-v ; ‘ . . V-. ‘ ^ • '-V • * 

0004 ■ N I N T R $ X X 'XXX'-l X -/aX-X / :i ' ’ } ' •" 

no-05 ■ nRdu$ . : • T :; .: V-- ' X-XXXXX^'XX:'.-' ■ •' ■ 

0006 V ■ N 101$ ' V XiXXX VVV- . , , — ; . ■'•■•.' 

0007 -'i NI02$ .''VV.V: V/ v V :V; ' .-s .■ 

ooio v nwdus X X ■■■'.•' - 

O011* ... SORT . X-rXX ■■■■X : : >>; 0^' v - >-■ -r .-V r-s"- > ^SV'-V^ ' 

0012 A T A N2 V ' • > '■ : " ■ ■ V; : ■ • ■ l’’/: XX : -‘X * :• ’■ X\ v7'- 

0013 SIN • "X'-XX ' 'W/ '-V V- ;W VVV'VV. 

.0014 'COS X. X OV ■X : K-X\ • 

0015 NSTOP5 r .ty : v XX% X^XXXX^- 'XXX ' - 


STORAGE ASSIGNMENT (BLOCK, type. RELATIVE LOCATION. NAME > • : ; 

0000 01.0020 lOOF 0000 h 01C022 101F ^ 000(1 010032 102F 000- 

0000 010037 107F 0001 000024 111G 0001 000033' 120G 0 00* 

0001 000114 155G 0001 V 000147 l?OG 0001 000162 177G 000*’ 

0000 R 010016 AMEAN 0000 R 00766© . ASPECH 0000 R 003000 ASPECY 000' 

DO 00 R 000000 DATAY 0000 R 010017 TM 0000 I 010014 I 000 j 

0000 1 010012 NS I G 0600 R 006000 PSPECH 0000 R 002000 PSPECY 00 0.; 



00100 1* 0************** ■E-***<HHHHHH>#******#****#*****#**«*#*S*****#*#«r 

00100 2* C DIMENSION VARIABLES. ' ; 7 

00101 4*. ' • DIMENSION DATAY< 2. 51 2). PSPECYC 512) , A SPEC Y < 51 2) , . ■ 

00101 5* : " 1 DATAH* 2. 512) .PSPECH (512) , ASPECH ( 512) , - - 

0C101 6 * XX.". 2 PTS(8) . •' ' ■ r ; - 

00101 7* C *>****»**»*»« *■ ~ **■&***««■***** ****»■&**»»***.&*****«•*.***#•& -»««««• 

00101 8* C " ASSIGN FOLLOWING -- 

00 101 9* C *»«■* * a ««>»-? »»»»««« *m> »*ee ********** * ******* **** « ****** 

00103 10* L NpTS = 8 .T",.'- 5 - V:';'' 

00104 11* : NPWR s .512 r V Vv,/ ■■'■■ 

00105 12* ' NS I G = 2 56 , : ’ v ; V- Vi V, ■■■ /■' X'l ■ /• 

00106 13* . CONST = 5. 84 68.: ; ; : • > : ' ; v . >- • 

00106 14* 0 ** ** &*»*■»* - ** » ****:. 

00106 . 15* C READ IN POINTS-' ' 

00106 16* .c* 5V * *»* ******* •’ >*******************»*****#*«-*^*iH»*-*?nHHHJ*** - *'*."f 

00107 17* .. READ (5.100) < PTS ( I ) » I =1 . NPTS ) 

00115 18* 100 FORMAT (4f 16. 3) 

001.16 19* sum = o. o. 

00117 ' 20* DO 1 I = 1 . NPTS • ’ ' - v* ;' ‘ 

00122 . 21* SUM = SUM + PTSC I ) ’ ! 

00123 22* 1 CONTINUE ^ v' 

00125' • 23* . AMEAN = SUM / fLO AT (NPTS ) /;■ I 

00126 24* WRITEC6.101) . • ;C 

00130 25* 101 FORMAT (’THE FOLLOWING POINTS - MEAN WERE READ; * ) * • < 

00131 : 26* . DO ? ! = l.NPTS . • 

00x34 27* PTS< I ) = F TS < I ) - AMEAN * ; 


0 01/35 

00141 
30142 

00142 
00142 
00142 
00144 
00147 

00150 
0 C 151 

00151 
0D151 
00151 
00153 
0 0154 
00157 
00160 
00161 
00162 
00.163 
00163 
00163 
00163 

00165 

00166 
00167 

00172 

00173 

00174 
00174 
00174 
00174 
00176, 
00 20.1 
00202 
Q0203 
00204 
00204 
0 0204 
00204 
00206 
00207 
00212 
00213 
0021.4 
00215 
00217 
00221 
00222 
00224 

00224 

00225 
00230 

00235 

00236 

00240 

00241 


'?■ 


28* 


Write (6 

•102) I.PTS(I) 

29* 

102 

FORMAT ( 1 

POINT C * , 13, • )= * ,F16.3) 

v30* 

2 

CONTINUE 

7 • 

31* 

C* * * ** * 


32* 

C 

GENERATE 

y POINT ARRAY. 

33* 

C * * V 


-:00 0C0000000000»«0 S a. S oa«00000000 r 00 00000000*0* 

34* 


•DO 3 1 - 

1,8 

35* 


DATAY( 1 , 

I ) = PTS(I) 

36* 


D A T A Y ( 2 > 

i ) = o.ooooo . : • 

37* 

3 . 

CONTINUE 


38* 

c**> 

00 00 *000 0* 

1! ? 0 0 0 0 0 0 0 0 0 0 0 0 0 9 <> * 0 * « 0 «»« 3. * * 0 0 0 0 £ 0 * * 0 » * * 0 0 * * *0. # 

39* 

c 

COMPUTE 

PPT. POWER SPECTRA AND ANGLE SPECTRA OF Y. 

40* 

C* * - 

0000000000 

000000000000 00 00 0*0000 0 000000000000 0 0*00** * * * 


41* CALL FOUfU<DATAY,!\IPWR,+l) 

42* DO 5 I = l.NPWft 

4 3* PSPECYO) = DATAYCl, 1 )**2 + DATAY ( 2 » I )**2 

44* /: PSPECY ( I ). = SORT ( PSPEC Y( I ) ) 

45* -- PSPECH (!) - PSPECY ( I ) 

46* ASPECY(I)- = ATAN2(DATAY(2»I )» DAT AY (1*1)) 

47* 5 CONTINUE .. ... 

48* C **-********** c r: **************************** *********** ****** 

49* C SET POINTS FOR X, . 

50* C*'-^**^*******"' 5 ********************************************* 

51* / AS PECK ( 1 ) = 0.00000 .7 

52* Fm = 0.00d00 - /; - ; / L - ■ 

53* : DO 6 I = 2 ,NPWR j - 7/-V: '//. 7, 

54* : ? :" " fm = fm + const ; •• ;v • 

55* 7/. ASPECH(I) s S I N ( FM ) ' 

*56* 6. CONTINUE -•/•/ .-.••• ' 

'• 57* c *® ******** ************************************** «■ 

58* C COMPUTE POINTS FOR H. 

59 * C* ° ***************** ** ** * * * ***** ****** ** ** ** ** ***«.***«*> *« ** *! 

60* DO 8 I = 1 , NPWR 

61* : ■•••; . ASPECH ( I ) = ASPECYd ) - ASPECH(I) 

62* DATAH(IfI) = PSPECH(I) * COS ( ASPECHC I) ) 

63* DAT AH (2,1) = PSPECH(I) * S I N ( ASPECH ( I ) > 

64* 8 CONTINUE 

65* C*^****®^******' 1 : 1 ********* 0 *********************************** 

66* ' C COMPUTE INVERSE FFT OF H . 

67* C************* : ***************************** **********<; ****** 

68* CALL FOUftM DAT AH ,NPWR» -1) 

69* DO 10 I = 1 , NPWR - . r . 

70* PSPECH < n = DAT AH (1 » I )**2 + D AT AH (2 » I ) **2 

: 71* ' ./ PSPECH(I) = SORT ( PSPECH ( I ) ) 

72*;:, ASPECH (I ) = AT AN2 ( DATAHC 2, I ) » DATAH( 1 # I } ) 

73* 10 continue . 

74* WR I TE ( 6 » 1 07 ) • 

75* 107 FORMAT ( 2RX » r SPECTRA OF INVERSE DATA H* ♦ ) 

76* • WRITE (-6 * 105 ) . / 7 ^ 

.77* 105 FORMAT <ltfl,lX,' FREQUENCY 1 » 2X , * POWER SPECTRA 1 ,2X , 

78* • 1' ANGLE SPECTRA • ) 

79* DO 11 I = 1* NS I G 

80* WR1 TE < 6,106 ) I , PSPECH ( I ) , ASPECH ( I ) . * 

81* 106 ’FORMAT ( IX » I12*2X,F14.3, 2X»F14.3) 

82* 11 CONTINUE • 

83* STOP : 

84* . . -END- • : • 


END OF COMPILATION: NO DIAGNOSTICS, 


frM.AP>IN .HAPf.ABS 

MAp 0023-04/04-13:47 -<,Q) 


frXOT- . AB$ . 

THE FOLLOWING POINTS - MEAN KERB READ: 



POI NT { 

1 ) = 

.450 

PO INT < 

2 ) = 

• .268 

POI NT ( 

3 ) = 

v . 120 

point ( 

4 ) = 

-.001 

PO I NT ( 

5 ) = 

. -.101 

point< 

6 ) = 

-.182 

pointc 

7) = 

V : -.249 

POINT ( 

8) = 

-.304 



SPECTRA OF INVERSE DATA H 
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83 
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24 

1.198 
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.86 

.346 

25 
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1.785 

87 

' .29 4 

- 26 

1.623 

1*781 'j 

: 88- 

.257 

27 

1.917 

1,77 7 

89 

.228 

28 

2.295 

1.774 

90 

.206 

29 

2.793 

' 1.771 , 

; 91 

.189 

' 30 

3.470 

1.769 

!< 92 - 

.177 

31 

4.428 

1.767 

■93 

.167 

32 

5.860 

1.766 

i .94 

.160 

33 

8.171: 

. 1 .764 

i 95 

.156 

34 

12.369 

1.764 

•j' 96 

.154 

35 

21.887' 

' 1.763 

97 

1 - .155 

36 

62.372 

1.761 

98- 

.15 9 

37 

44.452 
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! 99 

*166 
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64.118 
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39 
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